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ABSTRACT 

This  paper  presents  a  unified  approach  far  the  accurate  extraction  of 
specific  contact  resistivity  (  pt )  for  ohmic  contacts.  Using  2-D  simula¬ 
tions,  which  account  for  the  current  flow,  or  crowding  around  the  contact 
window,  we  have  analysed  the  resistance  data  obtained  from  the  Cross 
Bridge  Kelvin  Resistor,  the  Contact  End  Resistor,  and  the  Transmission 
Line  Tap  Resistor.  For  each  particular  structure,  a  universal  set  of  curves 
is  derived  that  allows  accurate  determination  of  given  the  geometry  of 
the  structure.  The  values  obtained  for  pc  are  independant  of  the  test 
structure  type,  its  geometry  and  the  contact  area  The  data  suggests  that 
in  the  past  researchers  have  overestimated  pc,,  and  that  contact  resistance 
will  not  limit  device  performance  even  with  submicron  design  rules. 

INTRODUCTION 

As  MOS  technology  scales  down  into  the  submicron  regime,  it  hat 
been  believed  that  the  series  resistances  contributed  by  the  source/drain 
contacts  increase  much  more  rapidly  than  the  other  resistance  components 
[1J.  This  is  because  the  contact  resistance's  scaling  approaches  X"J  for 
small  contacts  [2],  where  X  is  the  minimum  feature  size.  To  determine  if 
contact  resistance  is  a  limiting  factor  in  the  next  generation  ULSI,  it  is 
paramount  to  obtain  accurate  values  of  the  specific  contact  resistivity  pc, 
the  physical  parameter  that  governs  the  interfacial  contact  resistance 
between  the  contact  material  and  the  diffusion.  There  are  three  test 
structures  commonly  used  to  extract  pt:  the  Cross  Bridge  Kelvin  Resistor 
(CBKR)  (3),  the  Contact  End  Resistor  (CER)  [4],  and  the  Transmission 
Line  Tap  Resistor  (TLTR)  [5]  as  shown  in  Fig.  1.  In  all  of  these  struc¬ 
tures,  a  current  is  sourced  from  the  diffusion  level  up  into  the  metal  level 
via  the  contact  window  and  a  voltage  is  measured  between  the  two  levels 
using  two  other  terminals.  The  contact  resistance  for  each  structure  is 
simply  this  voltage  divided  by  the  source  current.  It  is  important  to  real¬ 
ize  that  each  device  measures  the  voltage  at  a  different  position  along  the 
contact,  hence  the  resistance  values  measured  are  different,  and  must  be 
clearly  defined  and  distinguished  from  one  another.  In  this  paper,  they 
are  referred  to  as  Rt  (Kelvin),  R,  (end),  and  Rf  (front)  for  the  CBKR, 
CER,  and  TLTR  respectively. 

There  are  two  problems  which  make  it  very  difficult  to  extract  pc 
accurately  from  the  contact  resistance  measurements,  using  the  present 
theoretical  1-D  equations  (given  in  Fig.  1).  The  first  problem  is  that 
when  these  three  structures  are  used  to  measure  similar  contacts,  the 
results  often  yield  conflicting  estimates  of  Pc-  The  second  problem  is  that 
when  contact  resistance  is  plotted  against  area,  a  sublinear  behavior  is 
observed  for  CBKR,  instead  of  the  expected  inverse  linear  behavior  [5-7). 
The  extracted  values  of  pc  not  only  appear  to  be  area  dependent,  but  also 
a  function  of  diffusion  sheet  resistance  R,  even  though  the  active  surface 
dopant  concentrations  are  the  same.  This  is  a  serious  problem  since  the 
variations  are  often  more  than  an  order  of  magnitude.  Previous  work  [8- 
12]  has  attributed  these  phenomena  to  two-dimensional  (2-D)  current 
Bow,  or  crowding,  in  the  diffusion-tap  area  around  the  contact  window. 
Those  analyses  partially  explain  these  phenomena  but  do  not  explain  how 
to  overcome  the  problem.  Hence,  there  is  a  need  for  deeper  understand¬ 
ing  and  unified  methodology.  In  this  paper,  a  unified  2-D  model  is 
presented  which  circumvents  these  shortcomings  and  provides  an 
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explanation  of  the  two  phenomena.  The  novel  extraction  method  derived  ' 
from  the  unified  2-D  model  provides  values  of  pc  which  are  accurate  and 
self-consistent  independent  of  geometry  and  test  structure  types. 

TWO-DIMENSIONAL  MODEL 

For  relatively  shallow  junctions,  a  diffusion  can  be  characterized  com¬ 
pletely  by  its  sheet  resistance.  We  shall  concentrate  here  on  semiconduc¬ 
tor  to  metal  contacts.  Since  the  metal  sheet  resistance  is  much  lower  than 
diffusion  sheet  resistance,  metal  is  considered  to  be  an  equipotential 
plane.  Therefore  the  current  Bow  in  the  diffusion-contact  system  can  be  i 
described  entirely  by  the  potential  in  the  diffusion  ’sheet*  which  is 
governed  by  the  Helmholtz  equation: 

#V  &V  V  ... 

U  +  =  (,) 

in  the  diffusion  area  directly  beneath  the  contact  and  by  the  Laplace 
equation: 

$*£-• 

elsewhere.  The  transfer  length  t,  is  defined  as  (  p*  /  R$  )*;  a  length  that 
characterizes  how  far  the  current  travels  in  the  diffusion  area  beneath  the  ; 
contact  before  passing  into  the  metal  level.  For  all  three  test  structures,  a 
ratio  between  a  measured  potential  (  V*  )  and  the  source  current  ( / ) 
gives  a  resistance  R*.  The  resistance  ratio  R*  /  Rt  can  be  expressed  as  a  i 
function  of  l, :  i 

f  o> 

4s*  I 

where  d\  is  ■  line  perpendicular  to  the  current  flow.  Note  that  the 
denominator  represents  the  total  current  Sowing  into  the  contact  window 
divided  by  the  sheet  resistance.  Using  numerical  techniques  [13),  we 
have  solved  equations  (1)  and  (2)  to  find  V(x,y,/,)  for  a  wide  ranging  set 
of  test  structure  geometries.  Then  R‘  /  R,  is  evaluated  via  the  solution 
of  (3).  By  comparing  the  experimentally  measured  resistance  ratio  to  the 
one  generated  by  the  2-D  simulations  using  /,  as  a  parameter,  a  unique 
value  of  I,  can  be  extracted  and  the  accurate  value  of  pc  obtained  for 
each  test  structure.  The  use  of  l,  allows  contact  resistivity  and  sheet 
resistance  to  be  lumped  together  as  a  single  parameter,  which  reduces  the 
simulation  task  to  one  independent  parameter  for  a  specified  test  structure 
geometry.  The  problem  is  further  reduced  by  considering  only  the  ratio 
of  measured  resistance  to  sheet  resistance. 

EXTRACTION  OF  pe  USING  2-D  MODEL 

In  Fig.  2,  the  ratio  of  the  Kelvin  contact  resistance  Rt  to  R,  is 
shown  as  a  function  of  contact  area  for  the  CBKR  structure.  The  overlap 
8,  which  is  the  difference  between  the  contact  size  and  diffusion  width  is 
maintained  at  2.5  pm.  The  diagram  shows  the  sublinear  behavior  in  the 
experimental  data  and  an  accurate  prediction  of  it  by  the  2-D  model. 
The  1-D  model  predicts  only  a  linear  behavior  (R4  =  pt  / 11)  because  it  is 
unable  to  account  for  current  crowding  effects.  The  diagram  shows  that 
there  is  a  large  deviation  from  the  1-D  model  as  I,  is  decreased.  This  is 
due  to  a  parasitic  component  added  to  the  ideal  value  due.lo. 
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current  crowding  effects.  This  is  more  pronounced  for  large  contact 
areas  and  small  values  of  fP  For  low  p,  and  high  Rr  serious  crowding 
occurs  and  2-D  simulations  are  required  to  accurately  model  the  behavior 
of  the  data  as  shown.  Good  agreement  between  simulations  and  meas¬ 
ured  values  implies  accurate  extraction  of  different  /,  for  the  metal  (W, 

A),  PtSi)  to  NVP*1  (B,  P.  As)  diffusions.  The  results  are  summarised  in 
Table  1. 

In  Fig.  3  the  2-D  model  is  compared  to  the  1-D  model  for  the  con¬ 
tact  end  resistor  where  S  is  2.5  pm.  The  2-D  model  shows  that  the  1-D 
model  severely  underestimates  this  resistance.  This  is  due  to  current 
flowing  around  the  overlap  area  and  into  the  contact  end  The  error  can 
be  as  high  as  scvetal  orders  of  magnitude  for  low  values  of  /,  [14].  This 
shows  that  when  the  1-D  equation  for  the  CER  is  used  to  extract  low  pc 
values,  the  true  pc  is  overestimated  by  at  least  an  order  of  magnitude. 

To  study  the  effect  of  the  current  flowing  in  the  overlap  region, 
CER  test  structures  were  fabricated  in  which  the  diffusion-tap  width  w 
was  varied  while  the  contact  size  is  kept  constant  at  3  pm.  The  data  and 
simulations  are  shown  in  Fig.4.  For  curve  1,  the  value  of  R,  is  indepen- 
dant  of  w,  implying  that  the  1-D  model  is  sufficient  for  this  high  value  of 
pc  But  curves  2  and  3  show  a  strong  dependance  on  w,  which  cannot  be 
accounted  by  the  1-D  model.  Once  more,  the  2-D  predictions  track  the 
measured  data  accurately  as  shown.  For  PtSi  to  P  diffusion  (curve  2),  the 
extracted  value  of  pe  is  3  fl  pm2,  which  agrees  with  the  CBKR  value 
obtained  earlier  on  the  same  wafer.  We  have  also  extracted  an  identical 
value  of  pt  with  the  TLTR  structure,  and  it  is  much  lower  than  previ¬ 
ously  reported  values. 

The  third  structure  to  be  examined,  the  TLTR,  measures  the  "front" 
resistance  indirectly.  The  "front"  .potential-  the  contact  potential  near  the 
leading  edge  of  the  contact  is  larger  than  the  potential  on  the  side  as 
sensed  by  the  CBKR  and  much  larger  than  the  contact  end  potential  in 
the  CER  case.  Therefore,  the  2-D  crowding  effect  has  much  less  impact 
on  the  TLTR  than  the  other  two  structures.  In  fact,  it  is  possible  to 
design  a  TLTR  which  is  essentially  1-D  by  making  the  contact  width 
much  larger  than  its  length.  The  effect  of  2-D  crowding  on  the  front 
resistance  of  square  contacts  is  shown  in  Fig.  3.  For  these  contacts,  the 
1-D  model  is  off  by  less  than  an  order  of  magnitude  even  for  very  low 
values  of  f„  signifying  that  the  current  crowding  effect  is  not  serious  for 
this  structure.  Unfomnately.  the  TLTR  structure  introduces  another 
source  of  error  because  it  can  only  measure  the  front  resistance 
indirectly.  The  measured  resistance  is  the  sum  of  the  front  resistance  and 
series  diffusion  sheet  resistance,  and  these  two  components  must  be  iso¬ 
lated  by  extrapolation  techniques.  i 

The  graphs  displayed  to  this  point  are  useful  for  extracting  accurate  j 
values  of  pe.  By  making  several  devices  with  different  geometries,  the  ! 
Re  I R,  ratios  can  be  plotted  on  the  appropriate  curve,  and  pc  can  be 
extracted.  The  consistency  of  the  extracted  pe  values  can  and  should  be 
used  as  a  gauge  to  judge  the  validity  of  the  experiment.  So  far,  only  one  ! 
geometric  variable  can  be  varied  at  a  time  in  each  study.  In  order  to 
make  possible  the  extraction  of  p,  for  a  very  wide  range  of  design  rules, 
this  constraint  is  removed  in  the  next  section.  I 


GENERALIZED  EXTRACTIONS 


If  the  dimensions  of  the  test  structures  are  normalized  by  dividing 
each  by  a  characteristic  length  L,  then  equations  (1)  and  (2)  will  become: 
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The  general  resistance  ratio  in  equation  (3)  will  still  retain  the  same  form 
but  the  independent  parameter  is  now  -j-: 
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ALL  MAT:  'UAL  IN  THIS 


This  shows  that  there  is  a  simple  scaling  law  which  may  be  used  to  nor¬ 
malize  all  of  the  above  simulations: 


R'  ., 
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where  f  and  w  are  the  contact  length  and  diffusion  width,  respectively. 
The  characteristic  length  L  can  be  any  one  of  the  geometric  variables 
used  to  define  the  dimensions  of  the  contact  structures.  A  convenient  L 
would  be  the  overlap  S  because  it  is  not  varied  over  a  large  range. 

This  scaling  law  allows  all  of  the  previous  simulations  to  be  nor¬ 
malized,  providing  three  sets  of  generalized  curves.  (It  is  still  necessary 
to  differendate  between  front,  end,  and  Kelvin  resistance.)  These  curves 
eliminate  the  need  for  further  2-D  simulations  since  all  practical  dimen¬ 
sions  of  contacts  and  diffusions,  and  all  contact-resistivites  and  diffusion 
sheet  resistances  are  contained  in  one  diagram.  The  application  is 
demonstrated  in  Fig.  6.  The  resistance  ratio  is  plotted  against  the  dirren- 
sionless  ratio  /  /  S  with  parameter  I,  I  8.  When  l, »  8,  the  behavior  is 
essentially  1-D  and  the  curves  are  straight  But  as  l,  decreases  below  6, 
the  2-D  current  crowding  effect  becomes  significant  especially  in  the 
cases  of  large  contacts.  The  extraction  accuaracy  in  this  case  is  rather 
low  as  is  shown  in  the  close  packing  of  the  curves  in  this  regime  of  low 
Pc-  As  contact  size  shrinks,  the  curves  spread  out  more  evenly  and  the 
extracdon  accuracy  is  improved.  Also  note  that  the  effect  of  reduced  8 
results  in  improved  accuracy. 

The  generality  of  the  universal  curves  is  demonstraicd  in  Fig.  7. 
The  2-D  simulations  and  experimental  data  of  two  sets  of  CBKR’s  fabri¬ 
cated  on  the  same  wafer  with  only  a  difference  in  6  (1.23  pm.  2.5  pm) 
are  shown.  The  extracted  /,  is  the  same  (0.3  pm),  indicating  once  more 
that  pe  is  independent  of  geometry  of  test  structure. 
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A  similar  set  of  universal  curves  for  CER  is  shown  in  Fig.  8.  The 
main  difference  is  that  the  values  of  the  contact  resistance  are  lower  than 
those  of  CBKR.  This  requires  measurement  of  significantly  smaller  vol¬ 
tages  if  current  densities  are  kept  the  same.  Furthermore,  the  resistance 
ratio  decreases  rapidly  as  8  increases  and  thus  the  sensitivity  of  the 
varaiadon  in  6  is  high.  This  is  because  the  end  contact  potendal  is  a 
strong  function  of  /,  /  /. 


The  TLTR  is  far  less  sensitive  to  8  because  it  detects  the  front  con¬ 
tact  potential.  Although  the  solid  lines  in  Fig.  9  indicate  that  the  TLTR 
has  slightly  less  geometric  sensitivities  than  the  CBKR.  there  is  another 
source  of  error  caused  by  the  indirect  measurement  method  in  extrapola¬ 
tion  of  the  contact  resistance  from  the  resistance  between  two  close  con¬ 
tacts  as  shown  in  Fig.  1.  When  R,  is  large,  the  extrapolation  becomes 
quite  difficult  and  sensitive  to  small  variations  in  the  separation  distance 
between  the  contacts.  But  the  next  example  will  demonstrate  that  with 
careful  electrical  and  optical  measurements  the  universal  extraction 
curves  yield  an  accurate  value  of  pc.  The  Rf  /  R,  values  of  the  TLTR 
devices  on  a  single  wafer  were  measured  using  separation  distances  and 
overlaps  measured  by  optical  and  electrical  techniques.  The  data  was 
then  plotted  as  symbols  along  with  the  universal  curves  in  Fig.  9.  At  first 
glance,  the  data  points  appear  scattered  and  poorly  fitted  to  any  single 
curve.  But  by  estimating  the  value  of  1,1  &  and  then  multiplying  this  by 
the  measured  value  of  8,  these  data  points  all  give  pc  values  within  20% 
of  each  other.  This  agreement  further  strengthens  the  validity  of  the  gen¬ 
eralized  method.  Using  this  technique  we  have  measured  pc  for  PtSi.  W 
and  A1  contacts  to  AT  and  P*  Si  for  a  wide  range  of  dopant  concentra¬ 
tion.  The  datailed  results  will  be  presented. 


CONCLUSIONS 


We  have  developed  a  technique  which  allows  accurate  and  unambi¬ 
guous  determination  of  contact  resistivity  even  at  low  values.  In  the 
examples  illustrated,  this  technique  repeatedly  yields  the  same  value  of  p , 
for  a  wide  variety  of  contact  dimensions  and  test  structures  on  a  single 
wafer,  fhe  results  of  the  2-D  model  are  presented  in  a  universal  form  so 
that  they  may  be  used  for  extractions  without  performing  any  additional 
simulauons.  The  new  model  points  out  that  the  I-D  model  seriously 
overesumates  the  specific  contact  resistance  due  to  2-D  current  crowding 
in  the  overlap  region.  This  implies  that  most  reported  values  of  pc 
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extracted  by  using  the  1-D  method  are  overestimated,  and  that  pc  will  not 
be  a  limiting  factor  for  ULSI. 
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Fig.  2  The  ratio  of  Kelvin  contact  resistance  to  diffusion  sheet  resistance  plotted 
as  a  function  of  contact  area  for  different  transfer  length  /,  for  CBKR  with  over¬ 
lap  being  2.5  jam.  Symbols  denote  measured  data,  solid  lines  are  simulations. 
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Contact  lenglh  ( i*m) 

Fig.  3  Simulated  end  contact  resistance  for  the  CER  structure  with  square  con¬ 
tact  holes.  Lines  and  symbols  both  denote  simulations  only.  Solid  liDes  arc  the 
2-D  model;  dashed  lines  are  the  l-D  model. 
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Fig.  1  The  three  test  struc¬ 
tures  used  for  measuring 
contact  resistance.  (a) 

Cross  Bridge  Kelvin  Resis¬ 
tor  (CBKR).  (b)  Contact 
End  Resistor  (CER).  (c) 
Trans  mission  Line  Tap 
Resistor  (TLTR).  The 

one-dimensional  (1-D)  for¬ 
mulas  for  extracting  specific 
contact  resistivity  are 

shown. 


s  _ ;n 

1  , 

VR.o,  ( 

■  ■  1 

w  tanh  (  i  p,  I )  j 

I  k 

ALL  MATERIAL  f'J  TMlS  3PAC 

[in 


4  *  21 J  pm 
Pc  =  1.1E-4  Clcm~ 


1,  =  6.4  pm 
pc  =  5  £-8  Clcmz 


li  =  4.3  pm 
Pc  *  6.9£-8  Qcm2 


— |-*-5  pm 


30  40 

W  ( pm) 


Fig.  4  End  resis¬ 
tance  as  a  function 
of  diffusion  width 
for  constant  con- 
lact  sire  oJ  5  pm. 
Symbols  denote 
actual  measure¬ 
ments;  lines  arc 
2-D  simulations. 
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Table  1  Summary  of  extraction  results  for  the  data  shown  in  Fig.  2. 
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Fig.  7  Illustration  of  the  use  of  the  CBKR  universal  curve  to  curact  /,  from 
structures  with  different  overlap  6  on  the  same  water  /,  =  0  5pm  for  ail  points. 
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Fig.  5  S  imulated  front  contact  resistance  for  the  TLTR  structure  with  square  con¬ 
tact  holes  Lines  and  symbols  both  denote  simulations  only.  Solid  hoes  are  the 
2-D  model;  dashed  lines  are  the  1-D  model. 


